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Abstract 

We perform variable separation in the quasi-potential systems of 
equations of the form q = —A~^\7k = —A~^'S/k , where A and A are 
Killing tensors, by embedding these systems into a bi-Hamiltonian 
chain and by calculating the corresponding Darboux-Nijenhuis coor- 
dinates on the symplectic leaves of one of the Hamiltonian structures 
of the system. We also present examples of the corresponding separa- 
tion coordinates in two and three dimensions. 
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1 Introduction 



In recent years, a new constructive separability theory, based on a bi-Hamiltonian 
property of integrable systems, was presented. In the frame of canonical co- 
ordinates the theory was developed in a series of papers (see also the 
review article [0]), while a general case was considered in [|] - [ |10|] . 

In this paper we apply the bi-Hamiltonian separability theory to solve 
Hamilton- Jacobi (HJ) equations for a class of quasi-potential systems, called 
bi-cofactor systems, systematically studied in recent papers |Tl|]-p!5[]. 

The structure of the paper is as follows. In Section 2 we briefly describe - 
the one-Casimir bi-Hamiltonian separability theory in case of bi-Hamiltonian 
systems with quadratic in momenta constants of motion and relate this the- 
ory with the classical Stackel theory. In Section 3 we present basic facts 
about quasi-potential systems and their special subclass called bi-cofactor 
systems. Then, in Section 4, an explicit form of transformation to sepa- 
rated coordinates is derived, proving the separability of all nondegenerated 
bi-cofactor systems. We sketch the proofs of more important results of this 
section independently of the general statements of Section 2 in order to make 
this article more self-contained, but we also point out the places where the 
general theory and our calculations meet. Finally, in Section 5, we illustrate 
the obtained formulas by some examples. 

2 Stackel separability of one-Casimir bi-Hamiltonian 
chains 

Let us consider a Poisson manifold Ai of dim M. = 2n + 1 equipped with a 
linear Poisson pencil H^ = Hi — ^Hq of maximal rank, i.e. a pair of Poisson 
operators (tensors) Hj : T*M. —>■ TM. each of rank 2n such that their linear 
combination Hi— ^Hq is itself a Poisson operator for any ^ G R (the operators 
Hq and Hi are then said to be compatible). 

Definition 1 A scalar function h : ^ Ti is Called a Casimir function 
(or a Casimir) of a Poisson operator H acting in M. ifllodh = 0. 

Here and in what follows d : C°°{Ai) T*A4 is the operator of external 
derivative (gradient operator) in A4 and the symbol o denotes the composi- 
tion of mappings. If - as it often happens - the Casimir of the pencil Hg 
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is a polynomial in ^ of order n 



h^ = h^C + hn-iC-' + - + ho (1) 

the by expanding the equation 11^ o dh^ = in powers of ^ and comparing the 
coefficients at equal powers we obtain the following bi-Hamiltonian chain: 

III o dho = 
Hi o dhi = Hqo dho 

■: (2) 

Hi o dhn = Ho o dhn-i 
= Ho o dhn 

From (Q) it follows that the functions hi are in involution with respect to 
both Poisson structures. If additionally all hi are functionally independent, 
then the chain is a Liouville integrable system. 

Let us now consider a set of coordinates {Aj,/ij}"^^ and a Casimir coor- 
dinate c' = hn on Ai so that {Aj,/ij}"^]^ are canonical with respect to IIq, so 
that 





/ 





no = 


-/ 







* 






We can try to linearize the system @ through a canonical transformation 
(A, /i) — *• (6, a) in the form bi = /ij = where W{X, a) is a generating 
function that solves the related Hamilton- Jacobi (HJ) equations 

dW 

^ri^^~Qx'> = ""r^ r = 0, ...,n. (3) 

In general, the HJ equations (Q) are nonlinear partial differential equations 
that are very difficult to solve. However, there are rare cases when one can 
find a solution of (^) in the separable form 

n 

W{X,a) = J2w^{Ka) (4) 

i=l 

that turns the HJ equations into a set of decoupled ordinary differential equa- 
tions that can be solved by quadratures. Such (A, fi) coordinates are called 
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separated coordinates. In the (a, b) coordinates the flow d/dtj associated with 
every Hamihonian hj is trivial 

dtj dtj dtj 

and the implicit form of the trajectories X{tj) is given by 

dW 

bi{X,a) = = Sijtj + const, i,j = 1, ...,n. (6) 



Theorem 2 A sufficient condition for (A, /i) to be separated coordinates 
for the bi-Hamiltonian chain is 



Hi{X,jj,c') = fi{Xi,iJ.i), i = l,...,n, 



(7) 



where 



Hi{X, /i, c) = c'X'l + hn-i{X, /i, c')Ar + ••• + ho{X, /x, c') 

and fi{Xi,^i) is some smooth function of a pair of canonically conjugate 
coordinates Xi, fii. 

We will briefly sketch the proof here. The condition (|^) can be presented 
in a matrix form 



/ Ai ■ ■ ■ 1 \ 



\n \n-l 
'^2 ^2 



n \n— 1 



A" A 



( 



J 



\ 



hn-l 



\ ho J 



( /i(Ai,/ii) \ 

/2(A2,/i2) 



or wh = f, (8) 



which is called a generalized Stdckel representation. Multiplying the HJ equa- 
tions (^) written in the matrix form 

h = a, with /i = (c, /lo)"^ and a = (c, a„_i, oq)^, (9) 

from the left by w one gets wh = wa, or according to the condition 
/ = wa i.e. 

dW 

/i(A„^) =cA^ + a„_iAr +- + ao, i = l,...,n (10) 
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which imphes the existence of the separated solution (^) for the system of 
HJ equations (0). In consequence, the system sphts into a decoupled set 
of ODE'S 



dWi 
dXi 



cA" + ttn-iX" + ... + ao, i = 1,. . . ,n. 



which concludes the proof. 

Let us now restrict to a special case, when fi{Xi, fii) is quadratic in mo- 
menta /ij 



(12) 



(for motivation of this choice see the end of this section) i.e. to the case 
when Ai = T*Q x R, where Q is some Riemannian (pseudo-Riemannian) 
manifold and T*Q its cotangent bundle. Then the condition (|]) can be put 
in the form 



\n— 1 \n—2 
^2 ^2 



\n— 1 \n—2 



K-2 /2(A2)/U2 + fl'2(A2) - A^C 



(13) 



This is a system of linear equations with respect to hi. It has the solution 

n 

/i, = J]Gf(A)/i^ + \/,(A) + cp,(A), ^ = 0,...,n-l (14) 

k=l 

where 



fc=i 



Aj = rij^ilAj — Aj), and Pi(A) are the polynomials of order n in A, defined 
by the relation 



(A-Ai)(A-A2)...(A-A„) = XIpiA^ 



(15) 
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(i.e. Viete polynomials). It can be shown that in the language of Rieman- 
nian geometry Gi are second order contravariant tensors on Q; the tensor 
Gn-i represents a contravariant metric tensor, i.e. the inverse of a standard 
covariant metric g on Q : Gn-i = g~^, while the remaining n — 1 ones are 
Killing tensors. From ( |l3l) it follows that the constants of motion for geodesic 
case have the form 



(16) 



k=l 



k=l 



where yjf = A" ^ / fi{\) is a classical Stackel matrix. It is worth mentioning 
here that the theory of separation of variables on Riemannian manifolds has 



been developed by Eisenhart [jT6| who used the concept of Killing tensors to 
analyze integrals of geodesic motion. The intrinsic formulation of Eisenhart 's 
theory has been presented in [jl^ where the concept of Killing web was in- 
troduced, while the application of the theory to bi-Hamiltonian systems has 
been presented in |l8l . 

From (^) it can be shown that in (A, /i) coordinates the operators IIq 
and Hi attain the form 





I 





no = 


-I 







* 






(17) 





A 


dhn-i/dfi 


ni = 


-A 


-dhn-i/dX 




* 






with the diagonal matrix A = diag{\i^,,,\n) and where the symbol * denotes 
the elements that make the matrices skew-symmetric. In fact there exists 
a whole family of separated coordinates (A', /i', c') which preserve the forms 
( [T^ ) and ( p!7D of hi and 11 j and which are related to the set (A, /i, c') by a 
canonical (gauge) transformation 



A^ = ±Ai, /i^ = ±/ij + §i{Xi 



l,...,n. 



(18) 



where d is an arbitrary smooth function. In the separated coordinates the 
Poisson pencil 11^ and the chain @ can be projected onto symplectic leaf Sc' 
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of Ho (dim S = 2n) producing a non-degenerated Poisson pencil on Sc' of the 
form Of: 



= 01 — ^60, where 
^0 = 



/ 
-/ 



01 



A 
-A 



(19) 



is a nondegenerate Poisson pencil on Sc'- Hence, Sc' is a Poisson-Nijenhuis 
manifold where the related Nijenhuis tensor A/" and its adjoint A/"* 



Af = 0io 0-^ 



A 
A 



AT* = 0-' o 0^ 



A 
A 



(20) 



are diagonal (here we use the name of Nijenhuis tensor for a second order 
tensor with a vanishing Nijenhuis torsion). This motivates the following 



definition ||T9| 



Definition 3 The coordinates (A, /i, c') in which Uq and Hi have the form 
(12) (or, equivalently , in which the tensor M = 0io0q^ has the diagonal form 



^^)), are called the Darboux-Nijenhuis (DN) coordinates. 



The tensors J\f and A/"* are not equal, since Af act on the space of vector 
fields while M* acts on the space of one- forms on A4. Notice that Pi{X) are 
coefficients of minimal polynomial of the Nijenhuis tensor 



(det(AJ - J\f)Y/^ = det(AJ - A) = ]^(A - A,) = J^P^^' 



1. (21) 



i=l 



i=0 



Let us now perform an arbitrary, not necessary canonical, c'-independent, 
nondegenerate coordinate transformation (A, p) {(l,p), preserving a quadratic 
dependence on momenta in all /;,/;, = 0, . . . , n — 1. It can be shown that 
after such transformation, i.e. in the variables {q,p, c') the operators IIo and 
III will have the following form 





-G(g) 





no = 


G^{q) F{q,p) 







* 






(22) 



Hi 



Giq) 
-G^{q) F{q,p) 



Uq o dh. 
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and the minimal polynomial of the Nijenhuis tensor becomes 



(Aet(il-M)r^^-^^m^, (23) 

det G 

Obviously, in a real situation we start from a given bi-Hamiltonian chain 
(0) in arbitrary coordinates {q,p,c'), derived by some method, and find the 
appropriate transformation to DN coordinates (A, /i, c'). In a typical situation 
one starts with the Hamiltonians hi that written in the coordinates (g, p, c') 
are quadratic in momenta p. If the transformation to (A,/i,c')— variables 
is such that A depends on q only, then the functions fi must have the form 
([T^ ) . This is precisely the situation one encounters in the case of bi-cofactor 
systems. In the next sections we will apply the ideas of this section in order 
to separate this class of systems. 



3 Basic facts about bi-cofactor systems 

Let us consider the systems of differential equations (in the fiat space R") of 
Newton form: 

q = M{q) (24) 

where q = (gi, . . . , g„)"^, M(g) = (Mi(g), . . . , Mn{q))'^ where denotes the 
transpose of a matrix, = qi(t), with t G R being an independent variable 
and where dots denote differentiation with respect to t, so that g'j = (Pqi/dt^ 
etc. A simple lemma below focuses our attention on those equations of type 
( p^ ) that posses a quadratic in q integral of motion. 

Lemma 4 Let 

n 

E{q,q) = q^A{q)q + k{q) = ^ qiAij{q)qj + k{q) 

where A = A^ is a non-degenerated, nxn symmetric matrix with q-dependent 

entries. Then E is an integral of motion for the system (that is E = 0) 
if and only if the following two equations hold: 

diAjk + djAki + dkAij = for all i,j,k = 1,... ,n (25) 
2AM + Vfc = (26) 
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Here and in what follows di = d/dqi and Vk{q) = {dik, . . . ,dnk). This 
lemma can easily be proved by calculating the derivative E = dE/dt. 

The equation ( pSj) implies that the matrix A is a Killing tensor. We will 
restrict ourselves to a class of solutions of (|25|) that have the form 



A = cofiG) (27) 

with 

G = aqq^ + (3q^ + qP^ + 7 (28) 

where cof means the cofactor matrix (so that cof{G)G = det(G') or, in 
case when G is invertible, cof{G) = det(G)G^^), a is a real constant, (3 = 
(/?!, . . . , Pn)'^ is a column vector of constants and where 7 is a symmetric 
n X n constant matrix. It is easy to show that for = 2 it is the general 
solution of (P^D. For higher n the general solution of (^) depends on n{n + 
l)^(n + 2)/12 parameters (its basis can be found e.g. in pO|) while (p7|) has 



only l + n + n(n + l)/2 parameters and is therefore far from being a general 
solution. It turns out, however, that this particular solution has interesting 
properties that make it worth studying. It originates in a natural way when 
one considers a broad class of Poisson pencils of type (0) (see [|12|)- It leads 
to the notion of bi-cofactor systems that admit many differential-algebraic 
properties, among them interesting recursion formulas that allow to generate 
non-trivial bi-cofactor systems from a trivial (geodesic) flow (see below). It 
has been recently generalized to the case of Riemannian manifolds [^. In 
this more general setting it can be demonstrated that the matrix G given by 
(|28| ) is a conformal Killing tensor, which leads to the orthogonal separability 
of geodesic Hamilton- Jacobi equation |T^ . 



The equation ( P6D implies that whenever det(y4) 7^ the force M can be 
written in the quasi-potential form M = —^A'^Vk, which generalizes the 
usual potential case and which reduces to the potential case when A = |/, 
where / stands for the identity matrix. Clearly, in case when our system (P^ 
has the second - functionally independent oi E - integral of motion of the 
form 

E{q,q) =fA{q)q + k{q) 

with an invertible matrix A, then it can be written in a quasi-potential form 
in two distinct ways. It motivates the following definition: 
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Definition 5 A system of equations 

q = M{q) = -lA'\q)Wk{q) = -^A'\q)Vk{q) (29) 
where A and A are two linearly independent matrices of the cofactor form 

)-m)-- 

A = cof{G), G = aqq^ + (3q^ + q(3^ + ^ 



A = cof{G), G = aqq'^ + Pq'^ + qp'^ + ^ 

and where k = k{q) and k = k{q) are two scalar functions, is called a bi- 
cof actor system. 



The bi-cofactor systems were first studied in |[TT|,[|T^ and in |T^, wliere 



they were called cofactor pair systems. In this article we will deal exactly 
with this type of systems: we will show how to perform separation of variables 
for quasi-potential systems of the bi-cofactor form. 

An important property of such systems is that they actually admit n 
constants of motion that are quadratic in q 



Theorem 6 (Hans Lundmark, /[T^/j // the Newton system ^4) has a bi- 
cofactor form (WW, then it has n integrals of motion of the form 



E,{q,q)=fAi{q)q + ki{q), z = 0, . . . , n - 1 (30) 

where the matrices Ai are defined as coefficients in the polynomial expansion 
of cofiG + S,G) with respect to the real parameter ^: 

n-l 

cof{G + ^G) = J2M' 

i=0 

with Aq = cof{G), An-i = cof{G) and where ko = k and = k. Con- 
sequently, in case when all matrices Ai are invertible, such system can be 
written in a quasi-potential form on n distinct ways: 

q = M{q) = -^A-'Vh, z = 0,...,n-l (31) 
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Remark 7 In the notation as above we have of course Eq = E and E^-i = 
E. 

Of course, for a given pair of matrices G and G not every function k will 
have a counterpart k that will satisfy the equation p9|). However, there 
exists a recursion formula that from a given bi-cofactor system produces a 
new bi-cofactor system. 

Proposition 8 Let q = M{q) be a bi-cofactor system of the form with 
the integrals given by Let also k^ = Y17=q ^j^* (with k^-i = k). Then 

the functions k, i = 0, . . . ,n — 1 defined as l^ = Y17=o through 

det(G + e6)- 
^ det(G) ^ ^ ^ 

satisfy the relation A'^Vk = Aj^Vlj for all i, j = 0, . . . ,n — 1 and are in 
consequence right hand sides of a new bi-cofactor system of the form I ^T\ ) 
(with k instead of ki). 



The proof of this statement can be found in |]T3|. This formula makes 
it possible to produce infinite sequences of bi-cofactor systems starting for 
example from a simple geodesic equation q = which is obviously of a bi- 
cofactor form with for example /cq = and = 1. Also, it can easily be 
inverted in order to express old quasi-potentials ki through the new quasi- 
potentials k: 

_l( det{G + ^G) \ 
^'--A det(G) ^"^^j ^^^^ 



In consequence, by using (|33|) it is possible to produce "lower" systems from 
the "higher" ones, and if we start from the geodesic equation g = 0, we obtain 
in general the family of "negative" systems different from the sequence of 
systems obtained from g = by the use of (^21). 

It has been proved [|1^], |T^ that the system ( pQ]) can be embedded in 



a bi-Hamiltonian system. In order to make this statement more precise, let 
us consider a following skew-symmetric operator pencil associated with the 
system (EOf): 
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Giq) 


P 


n5 = ni-eno = 


-G{q) F{q,p) 


M{q) + 2cN{q) 




* 






(34) 





-G{q) 







Giq) ~Fiq,p) 


-2ciV(g) 




* 






where both Hi and Ho are {2n + l) x [2n + l) matrices/operators acting in the 
(2n + l)-dimensional space Ai = R^""^^ with Cartesian coordinates labelled 

rp rp 

with c) where q = {qi-, ■ ■ ■ ■,qn) , P = {Pi, ■ ■ ■ ,Pn) , c G R. The n x n 
symmetric matrices G and G are exactly the matrices that define our system 
p9|). The n X 1 matrices and N are given by 

N = aq + (3, N = aq + i3. 

The n X n matrices F and F are defined by 

F = Np^ - pN'^, F = Np^ - pN^. 

As usual, the asterisk * denotes the elements that makes our matrices skew- 
symmetric. This operator pencil is a Poisson pencil precisely due to the 
fact, that the term M(g) can be represented as M(g) = — i^^^Vfcj for all 
i = 0, ...n — 1. One should also notice that this pencil is linear in the 
variable c and that it has a maximal rank 2n. Thus, according to Section 
2, the Casimir function (see Definition ^ obviously, in our variables d = 
{d/dq,d/dp,d/dc)^ ) of our Poisson operator 11^ is a polynomial of grad n 
in ^. In fact, the Casimir of 11^ has the form: 

h^{q,p,c) = p^cofiG + ^G)p + k^{q) - 2cdet{G + ^G) (35) 

where k^{q) = Y17=o Proposition || above. Thus, = 

Yli=o ^i^^ with the functions hi{q,p, c) given explicitly by 

hi{q,p,c) = Ei{q,p) - 2cA, i = 0, . . . ,n-l 

(36) 

hn{q,c) = -2cDn 



12 



with Di = Di{q) defined as 



5^A(g)e = det(G + eG) 



i=0 



so that Do = det(G) and Dn = det(G). The functions Ei in (|36[) are just 
the constants of motion (^) of our system (^). By expanding the equation 
o dh^ = in powers of ^ and comparing the coefficients at equal powers we 
obtain the hi-Hamiltonian chain of the form which by theorem of Magri 
is completely integrable in the sense of Liouville. 
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This means that the 



evolutionary equations 



d 
dt 



P 
c 



Hi o dhi 



Hn o dhi 



, n 



(37) 



associated with (H) are Liouville-integrable. 

Let us now investigate the relation of our chain (j^) (or of our Poisson 
pencil (p4D) with the system (^). One can show, that the last equation in 
(37) has at the hyperplane c = the following form: 



d 
di 



P 

c 



-2det(G') 



P 
M 





(38) 



d 


q 




p 


dt 






M 



so that the hyperplane c = is invariant with respect to this equation. On 
the other hand, if we set q = p in (^), we obtain its equivalent form 



(39) 



which differ from (|38| ) only by the coefficient — 2det(G). It means that both 
systems have the same trajectories in the (g,p)-space, although traversed at 
different speed. 

Lemma 9 (rescaling) Let x{t) = (xi(t), ... ,x„(t)) be a solution of a first 
order differential equation x = X{x) in R" with x(0) = Xq G R" and let 
a : R" ^ H be a scalar function. Then the function y{t) = x(r~^(t)) with 
r(t) = J -^^fjyjdt (chosen so that r(0) = 0) solves the differential equation 
y = a{y)X{y) with the initial condition y{0) = xq. 
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The proof of this lemma is elementary. This lemma implies that knowing 
a particular solution of (|38|) we can write down the corresponding solution 
of (p9|), just by identifying a with — 1/(2 det(G)). In the next section we will 
show how one can solve all the equations (0) (and thus also the bi-cofactor 
system (|39D ) by a procedure that separates variables in Hamilton- Jacobi 
equations that correspond to all Hamiltonians of the chain (^. 

More information about bi-cofactor systems as well as more detailed ex- 

In 



planations of the facts mentioned above can be found in [12|,|13],[15 
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and in one can find a generalization of many of the above mentioned 
statements to the case of Riemannian manifolds. 



4 Separation of variables 



In the preceding section we explained how our bi-cofactor system (p9D can be 
embedded in a bi-Hamiltonian chain (^. The main goal of this article is to 
present the procedure that leads to separation of variables in the Hamilton- 
Jacobi equations corresponding to the Hamiltonians hi. According to the 
rescaling lemma above, such a procedure will yield a solution of our bi- 
cofactor system as well. This procedure is based on the results obtained in 
and other papers. 



Let us begin by adjusting our coordinate system so that one of the new 
coordinates correspond to the foliation of M. into symplectic leaves of Hq. We 
will obtain it by rescaling c. Let us thus introduce the following curvilinear 
coordinate system in M.. 



Qi = Qi, Pi = Ph 



n, 



hn{q,c) = -2cdet{G) 



(observe that the hypersurfaces c = and c' = coincide). One should 
notice that this transformation does not depend on a particular choice of k 
and k (since the operator Hq does not depend on k and k) but only on the 
choice of G. In what follows we will write and Pi instead of q'^ and p'^. In 
the c'- variables the chain (H) has the same form; however, the explicit 
form of operators Hq and Hi changes to: 








-G 





Ho = 


G 


-F 







* 
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G 


-2pdet(G') 


ni = 


~G F 


2c' {N - GG-^N) - 2 det{G)M 




* 






(cf (0)), while the Hamiltonians hi attain the form: 



hi{q,p,c') = Ei{q,p) + d 
h„, = d 



I D. 



In the new variables the symplectic leaves of Ho have the desired form c' = 
const. Naturally, the question arises, if we could not choose a different pencil 
n so that the formula (BSp would not contain the factor — 2det(G), which 
would make Lemma |^ unnecessary, but we were not able to do so yet. 
We are now in position to present the main theorem of this article. 



Theorem 10 If the roots A, = Xi{q) of the equation 

det f G + Xg) = 



(40) 



are functionally independent, then in the variables A = (Ai, . . . A„,)"^, /i 
(/ii, . . ./i„)^; c' given by 



\i{q) as roots of 
^Q^cof (G + Xiiq)G')p 



^ = 1, 



n 



^n'^cof (G + Xi{q)Gj n 

where Q = GG^^N — N, the operators Uq and Hi attain the form 



(41) 



(42) 





/ 





no = 


-/ 







* 






(43) 





A 


dhn-i/dfi 


ni = 


-A 


-dhn-i/dX 




* 
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( cf (l^j ), with the diagonal matrix A = diag{Xi Xn), while the Hamiltonians 
hi have the form 

/.,(A,/x,c') = -ELitt^"%^ + ^'/^^W' ^ = 0,...,n (n>2) 

(44) 

hi{X,fi,c')=f{\,fi) + c'X, (n = 1) 

(cf ^1^)), with some functions fk depending only on one pair Xk, ^J'k of the 
variables, = nj^fel-^fc ~ ^j)' '^'^^ where Pi{X) are the Viete polynomials 
(fZ^/ Moreover, in the variables X,fi,c'the recursion formula attains the 
form 

1^{X) = deti^I - A)~k{X) - ^k^{X) (45) 

i.e. the formula ( ^) is invariant with respect to the change of variables 
q^X. 

Remark 11 It is worth mentioning that the coordinates X{q) defined by (^Tjj 
are in general not orthogonal, and that the gradients VAj(g) are eigenvectors 
of the matrix GG~^, i.e. (G + AjG')VAj = 0. Moreover, VAj(g) are G- 
orthogonal: VAjGVAj = for i 7^ j . In case when one of the matrices, 
say G, is equal to the identity matrix (that is when our system ( ^) becomes 
potential) the above transformation j^^-j^^) reduces to the classical formula 
for point transformation to separation coordinates for natural Hamiltonian 
systems. 



Remark 12 The terms Di/Dn in the Hamiltonians hi attain the form Pi{X) 
in precisely due to the fact that Di/Dn are coefficients in the polynomial 
expansion of det{G + AG) / det(G) = j^to Pi^' Tc/ 



Remark 13 The formulae (^^)-^^^) provide us with a transformation that 
is independent of the particular choice of the functions k and k in the bi- 
cof actor system (^^, i.e. this transformation will simultaneously separate 
all the bi-cofactor systems with the same matrices G and G. 



Theorem |T0| means that the coordinates (A, /i, c') are Darboux-Nijenhuis 
coordinates for our operators Ho, Hi. According to resuhs of Section 2, we 
have 
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Corollary 14 The Hamilton- Jacobi equations for the Hamiltonians hi{X, fi, d) 

= 0,1,. 



dX 



a.: 



where W{X, a) is a generating function for the transformation (A, /x) (6, a), 
separate under the ansatz W{X, a) = J2i=i ^ii^i^ ^'^^'^ system of ODE's of 
the form 



fk Xk 



dXk 



dXl + aiXl'^ + ■ ■ ■ + a„ 



(46) 



Proof. The Hamilton- Jacobi equations for Hamiltonians (|4^ ) can be 
treated as a system of n linear equations for functions /j. Applying the 
Cramer rule to this system we arrive at (^). ■ 

Thus, we are able to find W up to quadratures. If we denote the evolution 
parameter associated with hj by tj, then in the new variables (&, a) defined 
implicitly as usual: 



dW{X,a) 
dai 



dW{X,a) 

~~dxr~ 



1 n 



the fiow associated with hj has the trivial form 



so that the transformation 
(A, /i) (6, a) simultaneously trivializes Hamilton equations generated by all 
the Hamiltonians hi. 

We will now sketch the proof of Theorem Consider a symplectic leaf 
Sc' = {{lyPyC') '■ c' = const} of Hq. Let us choose a vector field transversal 
to the symplectic foliation of Ho as Z = It can be shown by direct 
calculation that 



LzUo = 0, LzUi = X AZ with X = Hq o d(Z(/i„_i)) 

where Lz is a Lie derivative operator in the direction of the vector field Z. 
The above relations guarantee that we can perform a projection of both Hq 
and Hi onto the symplectic leaf Sc' of Hq. The obtained 2n-dimensional 
Poisson operators 6o and 6i have the form: 







-G{q) 
-F{q,p) 



Giq) 
-G{q) F{q,p) 
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The corresponding Nijenhuis tensor Af = OioOq^ has the minimal polynomial 
of the form det(G'+ AG) / det(G) (cf (]23|)) and its roots are precisely the roots 
of pop. On the other hand, the roots of the minimal polynomial oi M define 
- according to (|21| ) - the first half of the transformation (g, p) — >• (A, /x) to DN 
coordinates in which the operators 9q and 9i have the form ([19D . Due to the 
last but one equation in the chain (|^) this implies that in the DN coordinates 
the operators IIo and IIi must have the form (0). 

We will now show that the remaining part of the transformation (g, p) 
{X, ^) to the DN coordinates, i.e. the expression for /i = fi{q,p), must be of 
the form (|i2|). We will do it in few steps. 



Lemma 15 In the DN coordinates the vector field X has a simple form: 

Proof. It is enough to calculate X = IIo o d{Z{hn-i)) in the DN- 
coordinates: 

d{Z{hn-l)) = d (^-^{hr._,)^ = d = -^^(Al + ■ ■ ■ + A„) 

where the last equality is due to the fact that Dn^i/Dn is precisely the same 
term that the last but one in the polynomial expansion of det(G'+AG) / det(G) 
that in turn is precisely the Viete polynomial p„_i = — (Ai + ■ ■ ■ + An). Thus, 

" d 

Ho o d{Z{K,_,)) = -Ho o d(Ai + . . . + A„) = 5^ — 



Lemma 16 The function 

n 
k=0 

(i.e. Casimir written in (A,/x,c') variables and evaluated at Xi) depends 



only on the i-th pair A,, fii of the variables A, /i i.e. 

Hi{X,fi,c) = fi{Xi,fii) 
for some function fi{Xi, fii). 
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Proof. We have that dHjdc' = ELo^idhk/dd = ELo^^^^^^/^n = 
det(G + AiG)/ det(G) = due to (|0]). For j ^ i we observe that dHi/dXj = 
fr + ^ifv + ^?fv + • • • + ^rlx"' z, J = 1, . . . , n. On the other hand, due 



to (M) and 



dhk dhn-i dhk dh 



so that 



dhk dhn-i I pk Pk-i Po 

I > o I ' ' ' I 



dX, dX, \X, X] xY' 



which substituted in the above expression for dHi/dXj yields zero. In a 
similar way one can prove that dHi / dpj = for i ^ j. ■ 

By Theorem ^ of Section 2 (and the pages that follow this theorem) this 
lemma means that our coordinates A, p indeed are separation coordinates for 



our systems. We will however continue our line of proof of Theorem 10. 



Corollary 17 From the above lemma it follows, by Cramer rule, that hi 
must have the form (j^/ 



Lemma 18 jl^/ Let X = li^d^Z [hn^i)) . Suppose that X^{Hi) = for some 



--2,3,.... and that X''{Hi) ^ for k = 1, . . . ,r - 1. Then 

p. = ■'^^"<^-> (47) 
In order to prove this lemma it is sufficient to integrate the relation 



X^'{Hi) = twice, using the fact that X = Y2^=i w~ ^^"^ Lemma |T6| and 



use gauge invariance of DN-coordinates in order to kill integration functions 
that appear after second integration. 

In our case the exponent r that "kills" Hi is equal to 3, since in our old 
coordinates {q,p, c') the vector field X has the form X = E - (...) ^ and 
since hi are quadratic in p. Explicit calculation of expression (|47|) for r = 3 
and in g, p, c'- variables yields exactly (|4^) . 

The recursion formula (^5]) is obtained by inserting / and A as —G and 



G in (1321). This concludes the proof of Theorem 10 
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5 Examples 



We will now illustrate the content of the presented theory with examples. It 
is worth to note that the general theory does not provide us with any tools 
for calculating the functions fi in (|44D (or in p6|)). Instead, we have to 
calculate these functions each time we perform the variable separation of a 
given bi-cofactor system. In case of bi-cofactor systems however it turns out 
that the functions fi always have the form (|l^). 

As a first example, let us consider the family of parabolic separable po- 
tentials introduced in [^. They have the form: 



h—n V / 



k=0 

with the other integral of motion given by 

= -q,qf + q,q,q, + qfv^^^^^ (48) 

We easily find the corresponding matrices G and G 



G 



-qi/2 
-qi/2 -q2 



G = \l (49) 



The recursion formula (p2D, applied to the geodesic equation g = with 
= ^ (i.e. with = 0, fci = 1) produces an infinite family of pairs of 
quasi-potentials k^\ k^(\ r = 0, 1, . . . (with k^^ = /cq, k^^ = ki) such that 
{cof{G))-^Vk^J^ = cof{G)Y^Vk^(^ = 2Vk'('\ The first few are: 

= -ql = -2q, 
ki^^ = 2q,ql kf^ = ql + m (50) 
kf^ = -qf - Aqlql A;f ^ = -Aq!q2 - 8ql 

where k^^ correspond to the potentials l^*-^^ up to a factor (—1)'' and where 
according to (|^) k^^^ = —q\k^(~^\ The corresponding potential-cofactor 
systems q = M^^\q) = —^{cof{G))^^Vk'^^ = —Vk^^^ have the form: 

q = M^^\q) = (0,2)^ 
g = M(2)(g) = (-2gi,-8g2)^ 
q = M(^\q) = (8gig2,4g? + 24gi)^ 
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so that the third one is already non-trivial. Similarly, by applying the formula 
(|33|), we can produce from the geodesic equation q = with = 1 (i.e. with 
k^^ = 1, kf^ = 0) the "negative" quasi-potentials: The first few of them are 
of the form 

kt'^ = 2q2/ql k[-'^ = -l/ql 

4"'^ = {^qI + qD/qI k[-'^ = -2g2/gt (51) 

kt^ = Aq,{2ql + ql)/ql k[''^ = -(4gi + q!)/qf 

and correspond to potential-cofactor systems g = M^^\q) = — | (00/(6*) )~"^Vfc| 
-Vk^{^ with 

q = M^-'\q) = {-2/ql0f 
q = M(-^\q) = i-8q,/ql2/qtr 

q = M^-^\q) = (-4(6g2 + qf)/ql 8q,/qff 

In order to check what variables will separate these systems, we have to solve 
the equation ( ^OD with G and G given as in (^91). An easy computation yields 

Ai(g) = q2- \lql + ql , ^2{q) = 92 + \fq[+~qi 

The formula (|4^) in this case reads 

V-Ai(g)A2(g) • 1 o 

Ai(g) 

and it is immediate to show that the above formulas present the classical 
point transformation to the parabolic coordinates. Thus, not only the poten- 
tials but even the corresponding chain (§) is separable in the parabolic 
coordinates which is perhaps what we should expect. After some algebraic 
manipulations, the above formulas can be inverted to 



qi = \/-AiA2 , g2 = - (Ai + A2) 
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_ A/-AiA2(/il - /i2) ^ _ jUlAi - /i2A2 

Ai — A2 Ai — A2 

which makes it possible to express the Hamiltonians hi in the DN coordinates 



(A,/i, c'). According to (p^) and ( pOD in the old variables {q,p,c') they have 
the form 

hi = li{q,p) + ki{q) + 

where 7j(q',p) = Ai{q)p is the geodesic part (geodesic Hamiltonian) of hi. 
Terms ^ = pi have in the DN coordinates the form of Viete polynomials 



(15) while the form of ki for a given bi-cofactor system obtained by recursion 



K^) can easily be established either by substituting the above expressions 



for gi(A) and Pi(A,/x) in (|50D-(pID or from the recursion relation (^51) . The 
result is 

(-3) _ -(Ai + A2)(A^ + A2) (-3) _ A^ + A1A2 + \l 

~ \3\3 5 ft^l - ,3,3 

(-2) _ -^1 + ■^1-^2 + A2 , (-2) _ — (Al + A2) 



\2\2 ' '^1 \2\2 



A1A2 A1A2 



A;?^ = AiA2 , A;i^) = -(Ai + A2) 
^ = -(Al + A2)AiA2 , fcf ^ = A? + A1A2 + Ai 
= AiA2(A? + A1A2 + Ai) , kf^ = -(A? + A?A2 + AiA^ + A^) 



The geodesic Hamiltonians 7^ have in our case the form 

^^2 1 , 2 I ^ 2 

7o = T ^ aAi^i + T r 2A2^2 

Al — A2 A2 — Al 



7l = T ^ ^Ai/ii + — iA2/i2 

Al — A2 A2 — Al 



22 



so that they indeed have the form with /j of the form ( ]T^ (as it has 
been pointed out at the end of section 2) and we can identify the functions 
fi{Xi) in (p!2| ) as fi{\i) = |Ai which is the form that can be used in order 
to solve the inverse Jacobi problem associated with equations (|46| ) and in 
consequence to separate the system (pOf ). 

As a second example we will consider a quite generic (but still two- 
dimensional) bi-cofactor system with matrices G and G of the form 



G 



(li(l2 ql 



G 



1 gi 
gi 2g2 



It is no longer potential. The first few of the bi-cofactor systems produced 
by the recursion formulae (|32D and (^) are defined by the quasipotentials 



(_2) _ 4 + 2g2 + 2g2 + q\ 



-2) _ -(2g2 -2g2 + g2^2 -4^2) 



4 



k 



(-1) _ 2 + g2 , (-1) _ 2g2 - 
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^(1) _ g2(2 + g2) 
2q2-q! ' ' 2g2-g? 



gi(2 + ga) ^(2) ^ gi(4 + 2^2 + qf + ql) 



and the corresponding forces M^^\q) are 

^^"'^ = A(?i(3 + g2),g2(4 + g2)f 
95 



(53) 



= ^ (gi,g2)' 



(g? - 2g2)^ 



(gi(l + g2),g2) 
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(g? - 2^2)^ 



{qi{ql + 2^2' + 4g2 + 4), q2{ql + 2^^ + 2^2))' 
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Also in this case the solutions of ( PH ) can easily be calculated: 



'™ 2(5^-292) '™ 2(gf-2,2) 

with A = 4g^ + (g2 — 2)^ > 0. The coordinate curves given by these equations 
consist of the non-confocal ellipses and hyperbolas ^6|] and an arbitrary point 



(9i;?2) in the g-plane may lay not only on intersection of an ellipse and a 
hyperbola (as it was the case in the classical separability theory) but also on 
intersection of two ellipses or two hyperbolas. The formulae (^2]) for fi{q,p) 
are in this case too complicated to be presented. The inverse relationships 
are however still quite compact: 

^V-AiA2(Ai + l)(A2 + 2) A1A2 



P2 



Ai + A2 + A1A2 Ai + A2 + A1A2 

^ _ ^ AiA2(/iiAi - //2A2 + /ii - /i2) 
Ai — A2 

-(Ai + 1)(A2 + 2) (/ii (A?A2 - A? + A1A2) - /i2 (AiAi - A^ + A1A2)) 



(Ai - A2)v/-AiA2(Ai + l)(A2 + 2) 
where we have chosen not to simplify the last expression since 
-AiA2(Ai + 1)(A2 + 2) = qUI/iql - 2g2)' 

is always non-negative. Applying these formulas we can express - after long 
algebraic manipulations - the geodesic Hamiltonians 7^ in the DN coordi- 
nates: 

70 = 4A?(Ai + 1)^? + 4A^(A2 + 

Ai — A2 A2 — Ai 



Ai — A2 A2 — Ai 

so that fi{\i) = 4Af(Aj + 1) in this case. The quasi-potentials (^) in the 



DN coordinates must attain the same form as the quasi-potentials (|50|)-(51) 
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do, namely the form given by ( |^ ) since the change of variables (^D|)-(^) is 
designed so that the recursion (|32|) in the DN coordinates always attains the 
form (|5l). 

In the end, let us consider a three-dimensional example with matrices G 
and G chosen as 





' 





1 " 






m2 


gi<?3 + 1 


G = qq^ + 





1 















1 










_ qm + 1 


52 gs 


ql 





' " 




" 





" 




" 





qi 


G = 





q^ + q[0 1 ] + 





1 










1 


q2 




1 













qi 


q2 


2q3 



The recursion formulas (^) and (|33D applied to the geodesic flow q = 



yield an infinite sequence of the quasipotentials k^^'^ = ^t=o^i^'^^^ ^ 

1 are 

i\ (7i ( 4- ^\ 

1 ^1 



^2 lW, 



. . . , -1, 0, 1, . . . .Some of them are 
, (-1) _ 2gig3 + 2gi + 1 _ gi(gi + 1) 



2gig3 + ^2 + 1 



2gig3 + ^2 + 1 



; "^2 



ql 



2gig3 + ^2 + 1 



(1) _ 2gig3 + ql + l (1) _ 2gig3 + 2gi + 1 (i) _ 2 + gi 
ftg - , fti - 2 ' 2 - ^ 

^(2) _ (2 + gi)(2gig3 + gi + l) , (2) _ 2 + gi(4g3 + 4 + 2gi - g|) 



gf 



ql 



, (2) _ gi - 2gig3 + 2gi + 3 

/Co n 



ql 



The corresponding forces M'^'^^ are 



(2gig3 + ql + 1)' 



, ?2(gi + 1) , qiqz - 1) 



M(^) = ^(0,0,1)' 

qi 
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= -1 {ql q^q^, q^q^ + gi + 3)' 



^^'^ = (liil^ + 4), qMqi + 3), qfqs + ql + 3gi + 6)^ 
Qi 

and they fast become complicated with the increasing \r\. In this case the 
formulas (^OD-(^) and their inverses are very complicated and can be handled 
only with the help of a computer algebra package. We will therefore quote 
here only the formulas for q{X) 

1 

qi = -2- 



q2 = 2 



Ai + As + A3 + 1 

v/-(AiA2A3 + A1A2 + A1A3 + A2A3 + Ai + A2 + A3 +T) 
Ai + A2 + A3 + 1 



?3 = 7 



1 Xi + Xl + Xj- 2(AiA2 + A1A3 + A2A3 + Ai + A2 + A3) - 3 
4 Ai + A2 + A3 + 1 



so that they are built of symmetric polynomials of order three in A. The 
quasi- 
form 



(r) 

quasi-potentials presented above attain in the above DN coordinates the 



^(-1) _ — (A1A2 + A1A3 + A2A3) _ Ai + A2 + A3 _ —1 

A1A2A3 ^ A1A2A3 ^ A1A2A3 



A;« = -A1A2A3 , = A1A2 + A1A3 + A2A3 , = -(Ai + A2 + A3) 



^0 ~ ('^l'^2 + A1A3 + A2A3)AiA2A3 

fcf ) = -(2A1A2A3 + AiA^ + AiA^ + X2XI + X2XI + A3A? + A3A^) 
= A^ + A2 + A3 + A1A2 + A1A3 + A2A3 
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which is in accordance with the recursion formula 
tonians 7^ have the following structure 



. The geodesic Hamil- 



7o 



A2A3 



4(1 + + ¥^ 4(1 + A^)/.^ ■ ■ ^ ^..^ 



A, 



A, 



4(1 + A3)/i^ 



7i 



%^ 4(1 + AO/.? + 4(1 + X,),l + Z%±M 4(1 + X,),l 



Ai 



A 



A, 



72 = ^ 4(1 + Ai)/x? + ^ 4(1 + A2)/i^ + + ^3)/^3 



so that they have exactly the form ( ^4|) with ( [l2l) and with /j(Aj) = 4(1 + Aj). 



One can see that in all the above examples the functions /j(Aj) do not 
depend on i, i.e. fi{Xi) = f{Xi). 



6 Conclusions 

In the present article we performed separation of variables for the recently 
discovered class of quasi-potential systems called bi-cofactor systems. These 
systems generalize the classical potential systems with additional, quadratic 
in momenta, integral of motion in the sense that they reduce to these systems 
in case when one of the matrices G and G is the identity matrix. In this spe- 
cial case the separation formulae (^D|)-(^2|) reduce to the well known form the 
classical separability theory formulas for separation of natural Hamiltonian 
systems by a point transformation. In the general case, however, these formu- 
las do not have the form of a point transformation and are to our knowledge 
new in literature. We concluded the article with some non-trivial examples 
in which the functions fi{Xi) in (|l^) actually do not depend on i and we can 
make a conjecture that it is always the case in bi-cofactor systems. 
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